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The electronic structures of Sri4_ :E Ca :E Cu2404i are cal- 
culated within the local density approximation. Around the 
Fermi energy there exist quasi-one-dimensional bands origi- 
nated from the ladder and chain layers. The nearest-neighbor 
inter-ladder hoppings are estimated to be 5-20% of the intra- 
ladder ones. Possible effects of Ca substitution on electronic 
structures and charge distribution are also discussed. 

PACS Number 71.15.La, 71.15.Mb, 74.25. Jb 



The doped spin-ladder compound, Sri4_ a; Ca a ;Cu2404i, 
has recently attracted wide attention as a possible new 
class of superconductor. Superconductivity_was observed 
at £=13.6 under high pressure 3.5~4 GPa.El The insulat- 
ing phase has a finite spin gap, and the superconductivity 
is expected to be driven by the spin-liquid ground state, 
as theoretically predicted in Ref. g. 

The superconductivity appears upon substituting iso- 
valent Ca for Sr, in addition to applying high pressure, 
whereas the formal Cu valence stays constant, 2.25 for 
all x. Considering this compound contains chain lay- 
ers as well as ladder layers, this indicates that the local 
charge distribution may change with Ca substitution and 
with pressure, and that superconductivity Hffds a special 
charge distribution. Various experimentscm and model 
calculations™] have shown that most of the holes exist 
on the chains for Sri4Cu2404i and that the Ca substitu- 
tion transfers the holes to ladder layers. In the present 
paper, we study the electronic structures and charge dis- 
tribution of Sri4_. E Ca : rCu2404i within the local density 
approximation (LDA).Q 

The crystal structures of Sri4Cri2404i and related 
compounds were reported in Refs. ^12]. They are com- 
posed of two subsystems, Sr 2 Cu 2 03 (ladders and Sr) 
and Cu02 (chains). r-Ihe chemical formula determined 
by structure analysialj is (Sr 2 Cu203)(Cu02); / with y = 
1.436. Due to interactions between the two subsystems, 
atomic positions of each subsystem are modulated by 
the periodicity of the other. Their structures were ana- 
lyzed by assuming a lacge unit cell with y approximated 
by a rati<mai,number,|3'EI or by the super space-group 
technique jlfJ'tHl If y is chosen as | or y, the chemical 
formula becomes SrioCui7029 and Sri4Cu2404i, respec- 
tively. 



We have calculated the electronic structures for 
M10CU17O29 and M14CU24O41 (M = Sr or Ca). As a 
starting point, the structure modulation is ignored for 
simplicity. Possible effects of the modulation will be 
briefly discussed later. The symmetry of M10CU17O29 
is chosen as face-centered orthorhombic F222 following 
Ref. ||. The a-axis is perpendicular to the ladders and 
chains inside layers, the 6-axis perpendicular to stacking 
layers, and the c-axis parallel to the chains and ladders. 
The lattice parameters are listed in Table I. c\ {c-i) is the 
fundamental period of the ladder (chain) part along the 
c-axis. 

Even with the rational approximation, the crystal 
structure of M10CU17O29 is still complicated and the 
unit cell contains 56 atoms. Such a complicated struc- 
ture is challenging for the ab initio calculation- We used 
the linear-muffin-tin-orbital (LMTO) methoofcj with the 
atomic-sphere approximation, since this is suitable for 
a large unit cell. In addition to 56 atomic spheres, 44 
empty spheres are inserted around Sr and Cu02 chains. 
The positions and size of the empty spheres are optimized 
by the method explained in Ref. [14| We have performed 
the self-consistent calculations with 63 k points in the ir- 
reducible Brillouin zone for M10CU17O29 and 21 k points 
for M14CU24O41. 

First, we show the total and partial density of states 
(DOS) in Fig. 1. The Fermi energy is set at Ep=0. 
The calculated DOS of Sri Cui7O29 and Sri 4 Cu2404i 
are similar to each other. This indicate that the small 
unit cell is enough to investigate the qualitative features. 
These compounds are calculated as a metal, leading to fi- 
nite DOS at Ep. To reproduce experimentally observed 
insulating behavior, it would be necessary to take into 
account electron correlations along with possible charge 
ordering. The states in the region — 7^2eV are mainly 
composed of the Cu d- and O p-orbitals. Both Cu and 
O partial DOS distribute broadly in this energy region, 
which indicates strong hybridization of these orbitals. 
The Sr s-orbitals slightly hybridize with the Cu d- and 
O p-orbitals, while their main weights appear at higher 
energy. 

The width of the valencej-band is consistent with the 
photoemission experiments.113 The two broad peaks ob- 
served in the photoemission spectrum at binding energy 
3eV and 5.5eV are actually seen in the present calcula- 
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tion. The former peak may be assigned to the peaks of 
the total DOS at — 2~— 3eV, and the latter peak around 
— 5eV. However, the position of Fermi energy is about 
leV lower in the calculation. This may be because the 
Fermi energy is not correctly calculated for this metallic 
ground state. 

The chain and ladder layers have different structures of 
partial DOS around Ep. The chain DOS forms a band 
with width of about leV, separated from the main va- 
lence bands with an energy gap of leV. On the other 
hand, the ladder DOS is smaller than the chain DOS and 
distributes in a wider energy region. The peaks at O.leV 
and 1.5eV in the ladder DOS are due to the edge sin- 
gularity of two quasi-one-dimensional bands, as we shall 
explain later. The inter-layer hybridization between lad- 
ders and chains seem to be small, since the peak positions 
of their partial DOS are not correlated. 

We show the energy bands near Ef in Fig. 2. The 
bands are most dispersive along the symmetry line T-Z, 
which is parallel to chains and ladders, and then along 
T-X (perpendicular to chains and ladders in layers). The 
small dispersion perpendicular to layers (r-Y) indicates 
that the interlayer hoppings are weak. We identify the 
character of each state by the weights of various atomic 
orbitals. Most of states near Ep have their weights only 
on either chain layers or ladder layers. This verifies small 
interlayer hoppings, consistent with the dispersionless 
character along r-Y. We also find that the chain bands 
near Ep are mainly composed of the antibonding combi- 
nation of Cu d^-orbitals and O p x - andp z -orbitals, while 
the ladder bands are composed of Cu d x i_ z 2 -orbitals and 
O p x - and p z -orbitals. 

We find that the calculated energy bands could be sim- 
ply understood in the following way. To see that, we plot 
in Fig. 3 the chain and ladder bands separately in the 
extended Brillouin zone. The states with more than 15% 
of their weights on Cu d x i_ z i -orbitals in the ladders are 
shown with large dots in Fig. 3 (a) and (b). The dif- 
ference is that those in (a) and (b) have an even and 
odd parity, respectively, concerning the mirror symme- 
try perpendicular to the a-axis. They both have the 
approximate periodicity of ^2L = izL ; corresponding to 
the fundamental ladder period. Thus, it is reasonable to 
interpret that these bands are generated by the folding 
of an energy band between k=(0,0,0) and (0,0,—). The 
band can be fitted by the following tight-binding disper- 
sion: 

£k= £o — 2<i cos(k z ci) — 2t 2 cos(2fc z ci) 
- [4if cos(ifc z ci) + At 2 cos(|fc z ci)] cos(|fc x a). (1) 

where t\ and t 2 are the nearest and second nearest neigh- 
bor hoppings between rungs in the ladders, while t± and 
t 2 between adjacent ladders. The results are summa- 
rized in Table II and the fitting is also shown in Fig. 3 
(a) and (b). With the five parameters for each band, we 
can reproduce the overall features. 

The inter-ladder hoppings are 5-20% of intra-ladder 



ones. These values suggest that the inter-ladder hoppings 
may not be negligible to discuss the electronic structures 
and physical properties. Experimentally, the anisotropy 
of resistivity p a / /x^is about 30 at temperature 100K for 
Sr3CanCu2404i,t£l larger than the present estimation, 
and it also has a large temperature dependence. There- 
fore it is natural to expect that the strong correlation 
effects such as hole pairing enhance the anisotropy, but 
we do not discuss it further in this paper. At high tem- 
peratures, the electron correlation may not be important 
and the prese-nt band anisotropy is consistent with the 
experiments.ll3 

The chain bands near Ep also have a pseudo one- 
dimensional character, while their dispersions are more 
complex as shown in Fig. 2(c). Here large dots repre- 
sent the states with more than 15% of their weights on 
Cu d xz -orbitals in the chain layers. We again fitted the 
dispersion by Eq. ([l]). The nearest neighbor hoppings, 
t\ , along the chains are smaller than those for the ladder 
bands. This is because of the near 90° angle of Cu-O- 
Cu bonds The second nearest neighbor hoppings, t 2 are 
largest and roughly twice the nearest neighbor ones. The 
sign and relative ratio of these two, t\ and t 2l are con- 
sistent with the semi quantitative estimation considering 
the spatial direction of atomic orbitals. 

Next, we discuss the effects of Ca substitution which 
corresponds to positive chemical pressure. The unit cell 
volume decreases and the lattice constant b, i.e., the dis- 
tance of subsequent layers, decreases most rapidly. We 
have calculated the electronic structure of fully substi- 
tuted Cai4Cu2404i with the lattice constants extrapo- 
lated from the available experimental data.El The inter- 
nal atomic positions are fixed to the same values with 
Sri4Cu2404i. The total and partial DOS have similar 
features to those of Sri4Cu2404i. The main change is a 
slight increase of the hybridization between Cu (i-orbitals 
and O p-orbitals. The enhancement of the hybridization 
is actually found as the wider total band width and the 
larger tight-binding hopping integrals shown in Table II. 
The increase is caused by the decrease of the atomic dis- 
tance. 

Let us now examine the hole distribution on the lad- 
der and chain layers, which is important to investigate 
the origin of superconductivity under high pressure. We 
use the occupation ratio of the chain bands near Ep as a 
measure of the hole concentration. The bands are com- 
posed of the antibonding combination of Cu d xz - and O 
p x - and p z -orbitals. Setting the O valence to —2, we de- 
termine that the valence of chain Cu would be +1 if these 
bands are fully occupied, and +3 if completely empty. 
This is exact if the interlayer hybridization were absent 
and is expected to hold qualitatively in the present case, 
too. 

With this assumption, we have calculated the chain 
Cu valence as P c = 3 — 2p, where p is the occupation 
ratio of the chain DOS to the whole number of states 
in this energy region. The valence of Cu on the ladder 
layers, Pe, is determined from the formal Cu valence P av 
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as Piri£+P c n c = P av (ni+n c ) where ni (n c ) is the number 
of Cu atoms on the ladder (chain) layers. The results 
are summarized in Table III. It is noticeable that the 
calculated Cu valence on the ladder layers is close to 2 
while that on the chain layers shows a large deviation 
from 2. It means that the most of the holes are on the 
chain layers. The Cu valence on the chain layers slightly 
decreases with the Ca substitution, indicating possible 
hole transfer from the chains to the ladders. However, 
the difference 0.02 ~ 0.05 is smaller than the proposed 
value, 0.18 for x = — > 11, from optical measurement .0 

So far, we have ignored structure modulations and pos- 
sible change of internal atomic positions upon Ca substi- 
tution. Since the largest effect on charge distribution is 
via Madelung potential, the detailed crystal structure is 
important to take into account. We have evaluated the 
effect of atomic displacements using the crystal structures 
where Sr atoms are uniformly shifted in the 6-axis direc- 
tion. If the shift is toward the chain layers, the potential 
for holes in the chain layersincreases and the holes would 
move to the ladder layers.! Our calculations show that 
the ladder Cu valence changes about 0.06 for = 0.13A 
shift of Sr atoms. 

Another important structure modulation is displace- 
ments of O atoms in the chain layers toward-the apical 
positions of Cu atoms in the ladder layers.EZl Since O 
atoms in the chain layers approach positively charged Sr 
and Cu, the potential for holes on the chain layers in- 
creases and the holes would transfer to the ladder layers. 
We have performed the calculation for Sr 10 Cui7O2g with 
the modulated structure reported in Ref. ||, with the lat- 
tice constants fixed. The ladder Cu valence increases to 
1.98, compared to 1.93 for the unmodulated structure. 
These results indicate the importance of structure mod- 
ulation to the hole doping in the ladders. 

In summary, we have studied the electronic structures 
of Sri4_2;CaxCu2404i within the local density approx- 
imation. The total density of states is compared with 
photoemission experiments with fairy good agreement. 
The energy bands near Fermi energy are from the ladder 
and chain layers. They have a pseudo one-dimensional 
character and the inter-ladder hoppings are 5-20% of the 
intra-ladder ones. The holes are mainly on the chains at 
x = 0. The change of the lattice constant is not enough 
to transfer the holes to the ladders. Detailed atomic dis- 
placements have important effects on the hole distribu- 
tion. 

We would like to thank to O. K. Andersen and 
O. Jepsen for providing us the TB-LMTO-46 program. 
A part of numerical calculations were performed on the 
NEC SX-3/4R at the Computer Center at the Institute 
for Molecular Science in Okazaki, Japan. 
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(a) Sr 10 Cu 17 O 29 



(b) Sr 14 Cu 24 41 
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FIG. 1. The total and partial density of states of (a) SrioCui7C>29 and (b) Sri4Cu24C>4i. 
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FIG. 2. The energy bands of Sr 10 Cu 17 O 2 9. r : (0,0,0), X:(2i ,0,0), Y:(0,^,0), Z:(0,0,2f ), and L:(|,f ,£). 
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(a) Ladder even (b) Ladder odd 

i i i I I i i i 




k=(0,0 I 0)-(0,0,27r/c 1 ) k=(0,0,0)-(0 I I 27r/c 1 ) 



(c) Chain 




k = (0,0,0)-(0,0,2^/c 2 ) 

FIG. 3. The energy bands of SrioCui7029 in the extended Brillouin zone. Large dots represent (a) states on ladder layers 
with even mirror symmetry, (b) those on ladder layers with odd mirror symmetry, and (c) those on chain layers. Solid lines 
are numerical fit. 
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TABLE I. Lattice constants (A) used for calculations. 
Sri Cui 7 O 2 4 Sri 4 Cu240 4 i Cai Cui 7 O24 Cai 4 Cu 2 404i 



a 


11.46 


11.46 


11.46 


11.18 


b 


13.40 


13.40 


12.40 


12.40 


c 


19.40 


27.65 


19.40 


27.10 


Cl 


3.88 


3.95 


3.88 


3.87 


C2 


2.77 


2.76 


2.77 


2.71 



TABLE II. Tight-binding parameters (eV) of the bands near Ef- The columns with ladder(+) and ladder( — ) are for the 
ladder bands with even and odd mirror symmetry, respectively. 





Sri Cui7C>24 




Ca 


10CU17O24 




Sri4Cu24041 




Cai4Cu24 04l 






ladder (+) 


ladder ( — ) 


chain 


ladder (+) 


ladder ( — ) 


chain 


ladder(+) 


ladder (—) 


chain 


ladder(+) 


ladder ( — ) 


chain 


£0 


-0.48 


0.35 


0.34 


-0.43 


0.43 


0.32 


-0.31 


0.46 


0.33 


-0.33 


0.54 


0.23 


ti 


0.42 


0.61 


-0.09 


0.44 


0.61 


-0.12 


0.41 


0.59 


-0.09 


0.44 


0.62 


-0.11 


t 2 


0.08 


0.08 


-0.15 


0.10 


0.08 


-0.20 


0.08 


0.07 


-0.17 


0.10 


0.07 


-0.20 


ti 


0.08 


0.03 


-0.03 


0.08 


0.04 


-0.04 


0.07 


0.03 


-0.03 


0.08 


0.07 


-0.04 


ti 


0.00 


-0.04 


0.03 


0.00 


-0.04 


0.04 


0.00 


-0.04 


0.03 


0.00 


-0.07 


0.05 



TABLE III. Estimated Cu valence. 



SrioCui7024 


Sri4Cu24C>41 


Cai()Cui7024 


Cai4Cu2404l 


average 2.235 


2.25 


2.235 


2.25 


chain 2.67 


2.64 


2.62 


2.62 


ladder 1.93 


1.97 


1.97 


1.99 
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